4 

Searching PAJ 



Kage l or z 



PATENT ABSTRACTS OF JAPAN 



(1 1 publication number : 2002-31 881 8 

(43)Date of publication of application : 31.10.2002 



1 — ■ 

(51)Int.CI. 




G06F 17/30 
606F 17/16 
G06T 7/00 




(21 Application number 
(22)Date of filing : 


2001-122755 
20.04.2001 


(71) Applicant 

(72) Inventor : 


CANON INC 

WASHISAWA TERUYOSHI 



(54) DATA PROCESSING DEVICE AND ITS METHOD, AND ITS PROGRAM 

(57)Abstract: 

PROBLEM TO BE SOLVED: To make data retrieval executable at a high 
speed by similarity based on an inner product relative to given vector data. 
SOLUTION: A sorting list relative to each component of a vector which is 
to be a prototype is made, and the nearest value to m-component values 
of test vectors X and -X are retrieved from the list, and respective 
positions are housed in PTR+ and PTR- (S303). A neighboring set to the 
number of k is initialized (S304), and a variable related to PTR is renewed 
(S306). Prototypes Yt are taken out successively, and determination 
processing is executed. If the similarity is larger than the minimum value 
up to the present, the prototype corresponding to the minimum value is 
replaced with Yt, and the neighboring set and the minimum value are 
renewed (S308-311). 
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* NOTICES * 
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l.This document has been translated by computer. So the translation may not reflect the original precisely. 
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3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 

[Claim 1] It is the data processor which extracts a prototype of a predetermined individual similar to the 1st given test 
vector from a prototype set described by multi-dimension vector as a retrieval result. A pretreatment means to create a 
list which arranged an identifier and a component value of each prototype of said prototype set in order of a component 
value about each component of said multi-dimension vector, and to calculate a square of a norm of each prototype, A 
test vector generation means to generate the 2nd test vector from said 1st test vector, A similarity count means to 
calculate similarity based on an absolute value of an inner product of the 1st test vector and each prototype concerned 
using a square of said list and a norm of each of said prototype, and the 1st and 2nd test vectors, A data processor 
characterized by having an abbreviation control means which omits future processings, and an updating means to 
update a retrieval result based on a value of said similarity, by decision of abbreviation conditions based on a value of 
said similarity. 

[Claim 2] A data processor according to claim 1 characterized by the above-mentioned test vector generation means 
generating said 1st test vector and a vector which a norm is equal and is reverse as said 2nd test vector. 
[Claim 3] When a dimension of a vector space as which said prototype and 1st test vector are expressed is set to d, said 
similarity count means If similarity is first calculated to 1 -dimensional subspace, said abbreviation control means 
judges that of abbreviation conditions based on similarity in the 1 -dimensional subspace concerned and these 
abbreviation conditions are not satisfied It is the data processor according to claim 1 which controls said similarity 
count means to raise a dimension of subspace and to perform similarity count, and is characterized by said updating 
means updating a retrieval result when not satisfying abbreviation conditions in addition, even if a dimension of 
subspace reaches d. 

[Claim 4] A data processor according to claim 1 with which said similarity count means is characterized by calculating 
similarity using the 2nd test vector and a thing which reduced a square of a difference of a prototype from a square of 
the 1 st test vector, a thing which reduced a square of a difference of a prototype, and a norm of a prototype from a 
square of a norm of a prototype. 

[Claim 5] A data processor according to claim 1 characterized by having a recognition means to recognize a class to 
which said 1st test vector belongs, based on a class to which a prototype of said predetermined individual extracted as a 
retrieval result belongs. 

[Claim 6] It is the data-processing method of extracting a prototype of a predetermined individual similar to the 1st 
given test vector from a prototype set described by multi-dimension vector as a retrieval result. A head end process 
which creates a list which arranged an identifier and a component value of each prototype of said prototype set in order 
of a component value about each component of said multi-dimension vector, and calculates a square of a norm of each 
prototype, A test vector generation production process which generates the 2nd test vector from said 1st test vector, A 
similarity count production process which calculates similarity based on an absolute value of an inner product of the 1st 
test vector and each prototype concerned using a square of said list and a norm of each of said prototype, and the 1st 
and 2nd test vectors, A data-processing method characterized by having an abbreviation control production process of 
omitting future processings, and an updating production process which updates a retrieval result based on a value of 
said similarity by decision of abbreviation conditions based on a value of said similarity. 

[Claim 7] It is the data-processing program which extracts a prototype of a predetermined individual similar to the 1st 
given test vector from a prototype set described by multi-dimension vector as a retrieval result. A head end process 
which creates a list which arranged an identifier and a component value of each prototype of said prototype set in order 
of a component value about each component of said multi-dimension vector, and calculates a square of a norm of each 
prototype, A test vector generation production process which generates the 2nd test vector from said 1st test vector, A 
similarity count production process which calculates similarity based on an absolute value of an inner product of the 1st 
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test vector and each prototype concerned using a square of said list and a norm of each of said prototype, and the 1st 
and 2nd test vectors, A data-processing program for making a computer perform an abbreviation control production 
process of omitting future processings, and an updating production process which updates a retrieval result based on a 
value of said similarity by decision of abbreviation conditions based on a value of said similarity. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[The technical field to which invention belongs] This invention relates to the equipment which searches or recognizes 
data by the similarity based on the inner product between the data expressed especially by the vector about the data 
processor which performs retrieval or recognition of data. 
[0002] 

[Description of the Prior Art] The inner product is widely used as similarity between the data expressed by the vector. 
For example, in character recognition or a voice recognition system, the sampled data is mapped to the characteristic 
quantity space stretched at the suitable base, and the data by which vectorial representation was carried out is 
memorized as a prototype. The inner product of the data and the prototype which were newly inputted is calculated, 
and input data is identified as a thing belonging to the class corresponding to the nearest prototype. Most, the bad count 
method of effectiveness is exhaustive search, and the computational complexity becomes the order of the dimension of 
a vector, and the product of the number of prototypes. 

[0003] it comes out in data base retrieval that the computational complexity of an inner product is recognized as a 
decisive failure. A data base can store now not only document data but an image and non-text data called voice by fast 
development of computer processing capacity in recent years. A keyword must be given in order to search non-[ these ] 
text data with the conventional keyword. If you want to avoid the time and effort of keyword grant, similar retrieval by 
the characteristic quantity vector must be performed. 

[0004] Moreover, even when searching document data, in order to enable more flexible retrieval, a document is 
expressed by the vector, and the algorithm which performs similar retrieval based on this vector is also realized. At this 
time, said computational complexity poses an essential problem of retrieval system implementation. The data number 
of cases stored in the usual data base exceeds hundreds of thousands of affairs. Therefore, it becomes the hopeless 
condition that computational complexity increases hundreds of thousands of times whenever one number of dimension 
of a vector increases. In order to avoid such a condition, the number of data which reduces the dimension of a vector or 
is calculated must be reduced. 

[0005] Since the former is equal to reducing the dimension of the presentation space of data, information required for 
data retrieval may not be enough expressed as a component of a vector. The latter is methodology which as for the 
number of data demanded as a retrieval result has a meaning when sufficiently small compared with the total of data. 
The problem dealt with by k-NN retrieval is this kind of thing, and some effective methods are proposed. 
[0006] k-NN retrieval is a method for searching for k near a test vector and identifying the class of a test vector based 
on those classes from the set of the prototype memorized by the system. In such a case, it is one of the important 
technical problems how k prototypes near a test vector can be found at a high speed. Such a demand exists also in data 
base retrieval. A retrieval user wants only some data near the search key which he specified among a lot of data stored 
in the data base, and no value is found out to other data, furthermore the value of an inner product, etc. The technology 
for meeting such a demand of a retrieval user is in agreement with the purpose of the high-speed algorithm of k-NN 
retrieval. 

[0007] When extracting k near a test vector from the set of a prototype, in order to mitigate the computational 
complexity of retrieval, it is common to structure the set of a prototype beforehand. Mitigation of the computational 
complexity of retrieval is expectable, so that it will carry out, if the property of data is made to reflect in the case of 
structuring. For example, in the case where a prototype is structured hierarchical, it is attained by repeating recursively 
actuation of dividing N dimension vector space where the prototype is expressed. 

[0008] It is the method of dividing that whose boundary used when carrying out division management of the space is a 
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hyperplane in K-D-B Tree [reference 1] and a rectangle field R-Tree It is SS-Tree about the method of dividing by 
[reference 2] and hypersphere. It is SR-Tree about the method of dividing space in the combination of [reference 3], a 
rectangle, and hypersphere. It is called [reference 4]. If N dimension vector space is furthermore rechanged into the 
space stretched by the characteristic vector about the covariance matrix of distribution of a prototype, more effective 
structuring is expectable to the computational complexity of retrieval [reference 5 and 6]. 

[0009] However, the application to the data with which the computational complexity and storage capacity for data- 
structure-izing increase exponentially to the increment in the number of dimension of a vector, and these methods are 
expressed by the high order former vector may be restricted as a matter of fact. Furthermore, to a fatal thing, the 
similarity based on an inner product is incalculable, if a test vector is not given. Therefore, when an inner product 
defines similarity, structuring of such data cannot be used. 

[0010] [reference 1] -- JT. Robinson: "The K-D-B Tree: A Search Structure forLarge Multidimensional Dynamic 
Indexes and" Proc. On ACM SIGMOD, pp. 10- 18 1981. [reference 2]A. Guttman: "R-trees: A dynamic index structure 
for spatial searching and" Proc. ACM SIGMOD, Boston USA pp.47-57, Jun.l984.[reference 3]DA. White andR.Jain: 
"Similarity indexing with the SS-tree," Proc. Of the 12th Int. Conf. On Data Engineering, NewOrleans USA pp.323- 
33 1 Feb [ 1996] [reference 4] Katayama, Sato: "Proposal of the index structure for the maximum contiguity retrieval 
to SR-Tree:high dimension point data" ****** (D-I), vol.l8-D-I no.8 pp.703-717, Aug. 1997. [reference 5] RF.Sproull: 
"Refinemants to Nearest Neighbor Searchingin k-Dimensional Trees and" Algorithmica 6 pp.579-589 991. [reference 6] 
DA. Hite and R. Jain: "Similarity Indexing: Algorithmsand Performance," Proc. On SPIE pp.62-73 1996. [001 1] On the 
other hand, there are also the loose structuring and the "loose" algorithm somewhat "wise" which has attained 
mitigation of computational complexity by the search algorithm which does not use a statistical property. This kind of 
algorithm is the only strategy which can be taken to the similar retrieval based on an inner product. Among those, one 
of the most fundamental things is the algorithm of Friedman and others called a projection algorithm [reference 7]. 
[Reference 7] Friedman F. Baskett and LJ. Shustek: "An Algorithm for Finding Nearest Neighbors" IEEE Structuring 
of the data demanded as pretreatment with JH. Trans. On Computers, pp.1000-1006, and a Oct. 1975. projection 
algorithm is sorting for every component of a vector, and is structuring based on a phase. That is, when a prototype is a 
d-dimensional vector, the sorting list of d pieces will be generated. 

[0012] Two kinds, the list Vj which stored j component value put in order by ascending order by this processing, and 
the list Ij which stored the corresponding prototype ID number, are acquired only for the number of dimension of a 
vector. That is, the n+lst values Vj (n+1) are beyond the n-th values Vj (n) from the head of Vj. Moreover, j 
component value Ylj (n) of the prototype Ylj (n) which makes Ij (n) an ID number, and (j) are in agreement with Vj 

fa) 

[0013] The principle of the projection algorithm for extracting 1 set of what is the closest to a test vector X from a. 
prototype set is explained using drawing 13 . Retrieval is performed along with the sorting lists Vm and Im of one 
piece chosen on suitable criteria. This supports having chosen m shaft among drawing. Since the data number sorting 
was carried out [ the data number ] by the component value is stored, as for Im, the list on a list is reflecting the phase 
on m shaft correctly. The value nearest to the m component X of a test vector X (m) is first looked for from Vm. This is 
set to Vm (j). The prototype corresponding to this is YIm (j). By a diagram, YIm (j) is Yl . Although Yl is the closest 
to X about m component, the distance in the whole space is not necessarily the nearest. 

[0014] Now next, the distance rho of X and Yl (X, Yl) is calculated. Then, only the prototype with which the value of 
m component belongs to an open interval (X (m) rho (X, Yl), X(m)+rho (X, Yl)) (the open interval A in drawing) may 
be closer to X than Yl, and it turns out that it is meaningful as an object of retrieval. In the example of drawing, the set 
of the prototype used as the candidate for retrieval is restricted further (X (m) rho (X, Y2), X(m)+rho (X, Y2)) (the 
open interval B in drawing) by inspecting the prototype Y2 near a degree about m component. Thus, it is a projection 
algorithm that mitigation of computational complexity is in drawing by making small the prototype set set as the object 
of retrieval based on the component value in 1 -dimensional subspace. 

[Problem(s) to be Solved by the Invention] Friedman's and others projection algorithm mentioned above is a method 
for the similar retrieval based on distance, and was not able to be applied for the similar retrieval based on the absolute 
value of an inner product, or a square. 
[0016] 

[Means for Solving the Problem] In order to solve a technical problem mentioned above, according to this invention, 
from a prototype set described by multi-dimension vector To a data processor extracted as a retrieval result, a prototype 
of a predetermined individual similar to the 1st given test vector A pretreatment means to create a list which arranged 
an identifier and a component value of each prototype of said prototype set in order of a component value about each 
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component of said multi-dimension vector, and to calculate a square of a norm of each prototype, A test vector 
generation means to generate the 2nd test vector from said 1st test vector, A similarity count means to calculate 
similarity based on an absolute value of an inner product of the 1st test vector and each prototype concerned using a 
square of said list and a norm of each of said prototype, and the 1st and 2nd test vectors, By decision of abbreviation 
conditions based on a value of said similarity, it has an abbreviation control means which omits future processings, and 
an updating means to update a retrieval result based on a value of said similarity. 

[0017] moreover, according to other modes of this invention, from a prototype set described by multi-dimension vector 
To a data-processing method extracted as a retrieval result, a prototype of a predetermined individual similar to the 1st 
given test vector A head end process which creates a list which arranged an identifier and a component value of each 
prototype of said prototype set in order of a component value about each component of said multi-dimension vector, 
and calculates a square of a norm of each prototype, A test vector generation production process which generates the 
2nd test vector from said 1st test vector, A similarity count production process which calculates similarity based on an 
absolute value of an inner product of the 1st test vector and each prototype concerned using a square of said list and a 
norm of each of said prototype, and the 1st and 2nd test vectors, By decision of abbreviation conditions based on a 
value of said similarity, it has an abbreviation control production process of omitting future processings, and an 
updating production process which updates a retrieval result based on a value of said similarity. 
[0018] By furthermore, the data-processing program which extracts a prototype of a predetermined individual which 
according to other modes of this invention gives and is similar to the 1st test vector of ******** from a prototype set 
described by multi-dimension vector as a retrieval result A head end process which creates a list which arranged an 
identifier and a component value of each prototype of said prototype set in order of a component value about each 
component of said multi-dimension vector, and calculates a square of a norm of each prototype, A test vector 
generation production process which generates the 2nd test vector from said 1 st test vector, A similarity count 
production process which calculates similarity based on an absolute value of an inner product of the 1st test vector and 
each prototype concerned using a square of said list and a norm of each of said prototype, and the 1 st and 2nd test 
vectors, A computer is made to perform an abbreviation control production process of omitting future processings, and 
an updating production process which updates a retrieval result based on a value of said similarity by decision of 
abbreviation conditions based on a value of said similarity. 
[0019] 

[Embodiment of the Invention] With the operation gestalt of** lst> book operation gestalt, the relational expression 
of an inner product and distance was drawn and the projection algorithm was used as high-speed calculus of an inner 
product. Furthermore, in order to evaluate an absolute value, the inner product to the given test vector and the inner ^ 
product to a test vector and the vector of hard flow were evaluated. Therefore, it is the hard structuring and the "hard" 
general-purpose thing which does not need the parameter which should be set up beforehand to a prototype set. 
[0020] Before giving details of the algorithm of this operation gestalt, the definition of a problem setup, a phrase, and a 

mark is clarified. . 
[0021] :omega= {Yl, Y2, YN} and Yj**Rd which set the set of the prototype Yj of N individual with which the 
problem to assume was expressed as a d-dimensional vector to omega ~ at this time, k large prototypes are extracted 
from omega about the absolute value of measuring rhoG (X, Yj) to test vector X** Rd given suitably. However, 
measuring rhoG (X, Yj) is : [0022] defined as an inner product. 



[External Character 1 ] 
p o(X,Yj) = XT G Yj = £ Z G< ™>»> X(m) Y i (n) (1) 

n-1 m=l 



Here, it wrote X (k) and Yj [ a test vector X and k component value of Prototype Yj ] (k), respectively. 

[0023] the feature of this operation gestalt is the about the same as measuring rhoG (X, Yj) ~ it is having adopted delta 

(Z, Yj) mentioned later as a function which gives a phase. It is equal to evaluating a prototype to descending of rhoG 

(X, Yj) and rhoG (- X, Yj) to extract a prototype to descending of the absolute value of rhoG (X, Yj). That is, if rhoG (- 

X,Vj) is large even if rhoG (X, Yj) is small, the absolute value of Yj of an inner product is large. 

[0024] • which can divide rhoG (X, Yj) and rhoG (- X, Yj) into two steps of following processings - Z=GXrhoG (X, 

Yj) =XTGY=(GX) TY=ZTY=rho (Z, Yj) rhoG(- X, Yj) -XTGY=-(GX) TY=-ZTY=rho (- Z, Yj) - rho (Z, Yj) and 

rho (- Z, Yj) are the inner products in orthonormal system here. On the other hand A ||2+||Yj||2-2rho (Z, Yj) top type is 

transformed further. :||Z-Yj||2=(Z-Yj)T(Z-Yj)=||Z which develops the square of the distance of Z and Yj and obtains a 

degree type - : which obtains a degree type - 2rho(Z, Yj)-||Z||2 - : which defines the right-hand side (or left part) of a 

=||Yj||2-||Z-Yj||2 top type as a new function delta (X, Yj) - now [ delta(- Z, Yj) =||Yj||2-||Z+Yj||2 ] like delta(Z, Yj) 
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=||Yj||2-||Z-Yj||2 Suppose that the absolute value of k prototypes and an inner product with those test vectors X is given. 
The minimum value is set to xi among the absolute values of these inner products. These k pieces presuppose that it is a 
recently contact in this time. 

[0025] Since Prototype Yj has the absolute value of an inner product smaller [ the absolute value of an inner product 
with the newly given prototypes Yj and X is smaller than xi or ] than all k prototypes if equal, or it is equal, it is 
guaranteed that it is not possible with k recently dots. That is, it can judge by whether a degree type is satisfied. : 
(xi>=delta (Z, Yj)) ** (xi>=delta (- Z, Yj)) A dot set can be updated (2), thus recently. On the other hand, :||Z-Yj||2 
>=||P(Z-Yj) ||2 in which the following inequalities are materialized about the 2nd term of the right-hand side of the 
definition type of delta (Z, Yj), however P are the orthogonal projection operators to subspace. Therefore, delta(Z, Yj) 
<=||Yj||2-||P(Z-Yj) ||2 (3) Similarly delta(- Z, Yj) <=||Yj||2-||P(Z+Yj) ||2 (4) A criteria (2) type is :(xi>=||Yj||2-||P(Z-Yj) 
||2) ** (xi>=||Yj||2-||P(Z+Yj) ||2) rewritten as follows using (3) and (4) types. (5) for example If the difference about m 
component of a vector, square ||Z(m)-Yj(m) \\2 of the sum, and ||Z(m)+Yj(m) ||2 are calculated and the criteria of (5) 
types are fulfilled, it turns out that Yj is not contained in a dot recently. 

[0026] (5) Evaluation by the formula is effective in respect of computational complexity rather than it uses (2) types. If 
the computational complexity for evaluating (2) types is actually estimated, since it is calculable beforehand, it is not 
necessary to take into consideration the 1st term of the right-hand side. As for the 2nd term, the multiplication for a 
number of dimension of a vector and addition are needed. On the other hand, the multiplication for a number of 
dimension of subspace and addition are enough for (5) evaluation of expression. Therefore, the more the dimension of 
the subspace where (5) types are materialized is small, the more mitigation of computational complexity is expectable. 
[0027] Hereafter, 1 operation gestalt of this invention is explained to details using a drawing. 
[0028] [Operation gestalt 1] The configuration of the data retrieval equipment of this operation gestalt is shown in 
drawing. 1 . I/O device 101 is equipped with a means to input retrieval data and commands, such as a communication 
line and a keyboard, and a means to output retrieval results, such as a communications control means and a display. For 
example, in the case of the computer of a stand-alone, the data which this I/O device 101 consisted of the keyboard and 
the indicating equipment, transmitted the input data inputted from the keyboard to inner product count equipment 102, 
and was transmitted from inner product count equipment 102 is outputted to an indicating equipment. 
[0029] On the other hand, in the case of the communication terminal connected to the communication line, this I/O 
device consists of a communication controller, and the input data inputted through the communication line is 
transmitted to similarity count equipment 102, and it transmits the data received from similarity count equipment 102 
to the address specified through the communication line. The set of the prototype of N individual expressed by the d- 
dimensional vector is stored in the data base 1 03 , and it is accessed by similarity count equipment 1 02. 
[0030] The procedure of the similarity count equipment 102>, next similarity count equipment 102 is explained using 
the flow chart of drawing 2 . 

[0031] First, it inspects whether there was any input from I/O device 101 at step S201, and if there is nothing, step 
S201 is repeated, and if it is, it will progress to step S202. At step S202, it inspects whether it is that to which input data 
updates the data of a data base, and if that is right, and that is not right, it will progress to step S203 to step S204. At 
step S203, pretreatment mentioned later is performed, and if it is completed, it will progress to step S201 . At step S204, 
it inspects whether it is the thing of computation, and if this input data is so, and it is not so, it will progress to step 
S205 to step S201 . At step S205, retrieval processing mentioned later is performed, and if it is completed, it will 
progress to step S201. 

[0032] Pretreatment performed at step S203 using drawing 9 is explained. 

[0033] With pretreatment, the sorting list about the square and each component value of a norm of the prototype 
belonging to a prototype set is created. Since the former only calculates and memorizes square ||Yj||2 of the norm of 
Prototype Yj, below, the latter processing is explained. A list of two kinds called the list Vj which stored j component 
value put in order by ascending order by this processing, and the list Ij which stored the corresponding prototype ID 
number is obtained only for the number of dimension of a vector. 

[0034] n is set as 1 at step S801 . At step S802, the pair of a n-component value and an ID number is created about each 
prototype of N individual, namely, ~ : - {(Yl (n), n) - (-- Y - two - (-- n --) ~ n --) ...} (YN (n), N) ~ a step - S - 
803 - the above-mentioned pair - the ascending order of a n-component value - rearranging . 
[0035] {(Yin (1), (n), In (1)) - (-- Yin (2), (n), and In (2)) - :Vn= which considers the list of a component value, and 
the list of an ID number as Lists Vn and In separately ... (Yin (N), (n), In (N)) and below} {Yin (1), (n), Yin (2) and 
(n), Yin (N) and (n)} (6) In= {In (1), In (2), In (N)} At the (7) step S804, n is incremented, if it is over the 
dimension d of a vector, processing will be ended, otherwise, it progresses to step S802. The relation of the list of two 
kinds is as follows. The n+lst values Vj (n+1) are beyond the n-th values Vj (n) from the head of Vj. Moreover, j 
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component value Ylj (n) of the prototype Ylj (n) which makes Ij (n) an ID number, and 0) are in agreement with Vj 
(n). 

[0036] The retrieval processing performed at step S205 using drawing 3 is explained shortly. As an input of retrieval 
processing, the number k of the prototype demanded as the vector X for retrieval (it is called a test vector below), and 
metric tensor G and a retrieval result is given. 

[0037] The index list of the base of vector space is created at the :Z=GX step S302 which multiplies a test vector by 
metric tensor from the left at step S301, and acquires Vector Z. :lambda= which this is a list which defines the sequence 
of the base which applies the rejection conditions mentioned later, for example, is a list corresponding to descending of 
the absolute value of the component value of a test vector X {lambdal, lambda2, lambdad} (8) -- : which writes L 
sets to be lambdaL from the one smaller again -- lambdaL= {lambdad-L +1, lambdad-L +2, lambdad} (9) - m is 
further set as lambda 1. 

[0038] At step S303, initialization processing of the variable PTR(ed) and related is performed. This processing is 
explained using dra wing 4 . 

[0039] The sorting list Vm of m component values is acquired at step S401 . Step S402 is searched for the value nearest 
to m component value X of a test vector (m) from Vm, and the location is stored in PTR+ at it. namely, : - |Vm(PTR+) 
X(m) |<=|Vm(j) X(m) | and **j** {1,2, N} -- similarly, it searches for the value nearest to -X (m) from Vm, and the 
location is stored in PTR-. namely, : -- |Vm(PTR-)+ -- at the X(m) |<=|Vm(j)+X(m) | and **j** {1,2, N} step 
S403 : which initializes a related variable as follows PTR+L=PTR+ 1, BND+L=0, and CAL+L=0PTR+H=PTR++ k 
near sets are initialized at 1, BND+H=0, CAL+H=0 PTR-L=PTR -1, BND-L=0, CAL-L=0 PTR-H=PTR-, BND-H=0, 
and the CAL-H=0 PTR=PTR+ step S304. This processing is explained using drawing 5 . 

[0040] At step S501, set NO (X) will be initialized to empty class soon. At step S502, t is set as 1 . At step S503, PTR 
later mentioned using drawing 6 is updated. Since the 1st term of the right-hand side of :second=Im(PTR) deltas=2|rho 
(Z, Yj) | which calculates the prototypes YIm (PTR) and delta s of ID number Im (PTR) at step S504, however a top 
type is calculated with pretreatment, only read-out from storage is required. 

[0041] step S505 - near set Nt- : which adds an ID number and the value of deltas to 1 (X) -- Nt(X) =Nt- 1 (X) + {(s, 
deltas) ~ if t is incremented at the} step S506 and k is exceeded, it will progress to step S507, otherwise, will return to 
step S503. 

[0042] At step S507, ID number s corresponding to xit-1 and it for the minimum value of deltas will be memorized as 
taut-1 in a set soon. Initialization of the near set to step S304 is completed above. 

[0043] At step S305, t is set as k+1. Renewal of PTR is performed at step S306. This processing is explained using 
drawing 6 and 7. 

[0044] At step S601, it is BND-L. 0 and PTR-L inspect whether it is 0, and if that is right, and that is not right, it will 
progress to step S602 to step S603. : which performs the following processings at step S602 ~ if it inspects, it comes 
out so and it is [ whether 0 and PTR-H have BND-H equal to PTR+H, and ] at BND-L=1 and the Dx-L=infinity step 
S603 - step S604 ~ otherwise, it progresses to step S605. : which performs the following processings at step S604 - if 
0 and PTR+L inspect how [ smaller than PTR-H ] it is, and BND+L comes out so and there is at BND-H = 1 and the 
Dx-H=infinity step S605 - step S606 ~ otherwise, it progresses to step S607. : which performs the following 
processings at step S606 - if 0 and PTR+H inspect whether they are two or more N, and BND+H comes out so and 
there is at BND-H=1, Dx-H=infinityBND+L=l, and the Dx+L=infinity step S607 ~ step S608 - otherwise, it 
progresses to step S609. : which performs the following processings at step S608 - at BND+H=1 and the 
Dx+H=infinity step S609, the product of four numbers, BND-L, BND-H and BND+L, and BND+H, inspects whether it 
is 1, if it comes out so and is, retrieval processing will be ended, otherwise, it progresses to step S610. 
[0045] At step S610, BND-L+CAL-L inspects whether it is 1, and if that is right, and that is not right, it will progress to 
step S61 1 to step S612. : which performs the following processings at step S61 1 - if BND-H+CAL-H inspects whether 
it is 1, comes out so and there is at Dx-L=|Vm(PTR-L)-Z(m) |CAL-L=1 step S612 - step S613 - otherwise, it 
progresses to step S614. : which performs the following processings at step S613 - if BND+L+C AL+L inspects 
whether it is 1, comes out so and there is at Dx-H=|Vm(PTR-H)-Z(m) |CAL-H=1 step S614 ~ step S615 -- otherwise, it 
progresses to step S616. : which performs the following processings at step S615 - if BND+H+C AL+H inspects 
whether it is 1, comes out so and there is at Dx+L=zVm(PTR+L)-Z(m) |CAL+L=1 step S616 - step S617 ~ otherwise, 
it progresses to step S618. : which performs the following processings at step S617 - if Dx-L is smaller than any of 
Dx-H, Dx+L, and Dx+H at Dx+H=|Vm(PTR+H)-Z(m) |CAL+H=1 step S618 - step S619 ~ otherwise, it progresses to 
stepS620. 

[0046] At step S619, the following processings are performed, and in :Dx=Dx-L, PTR=PTR-L, and CAL-L-0 step 
S620 which end renewal of PTR of step S306, if Dx-H is smaller than any of Dx-L, Dx+L, and Dx+H, and that is not 
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right, it will progress to step S621 to step S622. 

[0047] At step S621, the following processings are performed, and in :Dx=Dx-H, PTR=PTR-H, and CAL-H=0 step 
S622 which end renewal of PTR of step S306, if Dx+L is smaller than any of Dx-L, Dx-H, and Dx+H, and that is not 
right, it will progress to step S623 to step S624. 

[0048] The following processings are performed at step S623, and at :Dx=Dx+L, PTR=PTR+L, and CAL+L=0 step 
S624 which end renewal of PTR of step S306, the following processings are performed, and by update process 
of :Dx=Dx+H and PTR=PTR+H which end renewal of PTR of step S306, and CAL+H=0PTR, if the value of the 
variable relevant to PTR is changed and conditions are satisfied, the retrieval processing itself shown in drawing, 3 will 
be ended. 

[0049] By step S307, a degree type estimates (5) types by judgment processing later mentioned per drawing 8 at 
the :rho+=||Ys||-|| Z||rho-=||Ys||-||Z|| step S308 which calculates rho+ and rho-. 

[0050] At step S309, the return value of the processing performed at step S308 judges whether it is FALSE, and if it is 
FALSE, and that is not right, it will progress to step S3 10 to step S3 12. 

[0051] At step S3 10, the element corresponding to xit-1 will be deleted from a set soon. : which adds the prototype 

under current processing - Nt(X) <-Nt-l(X) {(taut-1, xit-l)}+ -- {(Im (PTR), Dx) -- at the} step S31 1, the minimum 

value of Dx of the elements of Nt (X) is stored in xit, the ID number is stored in taut, and it progresses to step S313. 

[0052] At step S312, a degree type is performed, t is incremented at t-:xi t=xi t-ltau t=tau 1 step S313 which progresses 

to step S3 13, if the result exceeds N, processing will be ended, otherwise, it progresses to step S306. 

[0053] The function count performed by judgment processing of step S308 is explained using drawing 8 . 

[0054] At step S701, j is set as 1. A degree type is judged at 2Dx <-Dx-(Ys(lambdaj)-X (lambdaj)) 2 step S703. : which 

performs the following processings at step S702 -- rho+=rho+-(Ys(lambdaj)-X (lambdaj ))2rho-=rho -- (Ys(lambdaj)+X 

(lambdaj)) -- If satisfied, TRUE will be set as a return value and processing will be ended, otherwise, it progresses to 

step S704. At step S704, if j is incremented and the result exceeds the dimension d of a vector, FALSE will be set as a 

return value, and processing will be ended, otherwise, it will return to step S702. 

[0055] Nt (X) when ending at step S3 1 3 is outputted as a retrieval result. 

[0056] The effect by the operation gestalt explained above is verified by computer simulation. 

[0057] [computer simulation] -- the number of prototypes demanded as a retrieval result in order to verify the 

effectiveness of the operation gestalt mentioned above -- k= 10 pieces and several prototypes -- computer simulation of 

retrieval was conducted to N= 10000 pieces, dimension [ of the : and the vector whose following items are experiment 

parameters ]: -- d= -- {-- : and CPU:PentiumIII whose item of the computer used for the 10 unit} experiment from 10 to 

100 is as follows, and - (500MHz) main memory: 128 MB-OS:Linux -2.0.36 -- in addition, programming language used 

C. 

[0058] [Experiment procedure] 

(1) Generate the prototype set which consists of a d-dimensional vector N individual using a uniform random number. 

(2) Use a uniform random number and it is d dimensions. One metric tensor of vector space is generated. 

(3) Generate one d-dimensional test vector using a uniform random number. 

(4) Perform exhaustive search. 

(5) Perform retrieval by the proposal algorithm. 

[0059] The average of the relative CPU time which shows the five above-mentioned procedures to a 10 times repeat 
and the following was calculated. 

Relative CPU time = (CPU time of a proposal algorithm) a /(CPU time of exhaustive search) experimental result is 
shown in drawingJO . The horizontal axis was shown for the dimension (number of dimensions) of a vector among 
drawing, and the axis of ordinate and the number N of prototypes were shown for the relative CPU time (CPU time 
ratio) as a parameter. 

[0060] Drawing shows that it is not concerned with the number N of prototypes, but the relative CPU time is increasing 
to the primary order with the increment in the dimension of a vector. However, the inclination at the time of N= 10000 
pieces is very small. The value of the relative CPU time at this time indicates the sufficiently small value to be- 
dimensional [ 1 00 ] or 1 8% also by high order origin 2% by ten dimensions. 

[0061] <Operation gestalt 2> The case where the similarity count mentioned above is applied to recognition equipment 
is explained. 

[0062] Drawing \ \ is drawing showing the configuration of the recognition equipment of this operation gestalt. the set 
of the prototype of N individual expressed by the d-dimensional vector is stored in the data base 103, and two or more 
prototypes exist to each category used as the candidate for recognition, respectively. Others are the same as that of 
drawing 1 . 
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[0063] In for example, face recognition, the d-dimensional vector stored in a data base 103 can constitute by connecting 
the person's information with the image in two or more angles of an everybody object (Mr. A, Mr. B, ...), expression, 
etc. In such a case, by processing described below, an input image and k near images are chosen, the frequency of the 
person number connected with each can be totaled, and the highest person of frequency can be specified as he is a 
person corresponding to an image. 

[0064] Moreover, the procedure of similarity count equipment 102 follows drawin g 2 like the operation gestalt 1. 
However, in step S204-205, retrieval processing turns into recognition processing which used retrieval. This 
recognition processing is explained using the flow chart of drawin g. 12 . 

[0065] In drawing .12 , step S301-313 are the same as that of drawing 3 . At step S314, the frequency of the class 
number related with the ID number which will belong to a set soon is calculated. An ID number is rearranged into 
descending of the frequency of a class number at step S3 15. Thereby, a class candidate is obtained by order with a high 
possibility that a test vector belongs. Among these, a class with the largest frequency may be outputted as a recognition 
result, and a predetermined individual or all classes may be outputted as a recognition candidate in order of frequency. 
[0066] In addition, even if it applies this invention to the equipment which consists of a single device, it may be applied 
to the system which consists of two or more devices. Moreover, the storage which memorized the program code of the 
software which realizes the function of the operation gestalt mentioned above may be supplied to equipment or a 
system, and you may attain by reading and performing the program code with which the computer in equipment or a 
system was stored in the storage. 

[0067] Furthermore, by reading and performing the program code with which the computer in equipment or a system 
was stored in the storage, also when an above-mentioned function is realized, it is contained by processing of OS which 
it not only realizes directly the function of the operation gestalt mentioned above, but is working on a computer based 
on directions of the program code. 

[0068] The storage which memorized the program code in these cases will constitute this invention. 
[0069] 

[Effect of the Invention] Since the similarity based on the absolute value of the inner product of a vector is calculable at 
a high speed according to this invention as explained above, based on this similarity, it is effective in the ability to 
perform retrieval of similar data and recognition of a class which belongs at a high speed to the given data by which 
vectorial representation was carried out. 
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40 itr, m&r^yz&zmmbtc. ek. «*mb* 

ff«r4fc©{c. ^^.6n/crx h^? h;Kc3*-T4rt 
Si. h^iit^lS]©^^ h^tc>Ft-r4rtS 

iSfflfffiOyt. {tit. XPhW:/|g£K:*t-r4" K 

yiffl6 t J^4 ) ©-r*4o 

[0020] *-mnm<D7>id y XA©i»^ai^4 
ijrtc. PpIsiss. ii^iia#©sa^Bjp>*HcL-cfc 

[0021] ®m?z>mm&> a-xyt^ h *t t-csi 

50 iSh/cNi©7'0 h $^7'Y,©*^4Qi-r4 : 



. • . (5) ftP32 00 2-3 1 8 8 1 8 

7 . 8 

Q = {Yi ,Y, • • ,Y» } ,Y, £ R d *T!£fl£ft£: 

C©£S> jl^£^x.<E>ftfcrX \"<f b)lXE R'tC. [002 2] 

ffffiPcCX.Y^ISttffiKfWbTA-^kffl©^ h£W mi) 

n=d m=d 

j oa(X,Yj) = XTGYj= £ J] G(ra,n) X(m) Yj(n) (l) 



-en-etxx(k). YtCiotniBbfc. 

[0 02 3 ] 439BftR©ttSttt. ftfiPc (X,Yi)i|lI 10 

3U7£«|ffi-f 5 Citt. P« (X,Y,)£p e (-X. 

t<>. Pc (X,Y 1 )*^?<"C*)Pc (-x,y, 

[0 024] p t (X,Y, )Rtf p« (-X.Y, )fi«T©J; 5 &2 
Z=GX 

P c (X,Y, ) = X T CY= (GX) T Y = Z 1 Y= P (Z , V, ) 20 
Po (-X.Y, ) = -X T CY= -(GX) T Y= -Z T Y= p (-Z.Y, ) 
C C V P (Z,Y, p (-Z , Y, )» iE«K SJRKte » 4F«J 

area*, -u. ztv x t<Dwmo2m*m.mi,'C'Jz£< 
*m : 

I Z-Y, 1 2 = (Z-Y, ) T (Z-Y )= I 2 I * + I Y, I * -2p (Z,Y, )* 

US$(Z,Y,)) A (|^S(-Z,Y,)) 



2p(Z,Y,)-IZ0 2 = |Y, I'- II Z-Y, I* 
±5$©:& M&Mt&mZffictm®. 8 (X,Y, )£ UT^ 

mz> : 

S(Z,Y,) = I Y, | '-I Z-Y, II 2 

S(-Z,Y,) = lY, I'-IZ+Y, I 1 
4. kffl©7"P h^-f^i. -ett6©rX F-^? F^Xi 

©^^©^fii^x.enTi^iTSo cn^rta© 

**HI©5 5. C©kffl#. $r$£ 

[0 02 5] Sr/c{c^e,ti/c7'P F *<f'7Yi ixi© 
-/7"Y,«. k(I©:£T©7'CI FW^JC*)^ F«9S©^ 

SiflWJ l^90t»©r. kH©*jSflfi&TW»)ff& 



C2) 



P(Z-Y, ) 1 



(3) 



5(Z,Y,)©^il^©*j2m2^{cB8Lr«T 
©*3£jfljBM-* : 

I Z-Y, 1 1 k J P(Z-Y, ) 1 1 30 

M£b. P»8U»SIBI^©itS3ti^fflSR-c*S. «£o 
*c * 

|Y, B 1 _ | PCZ-Y, ) 0 1 ) A 
WU*. hJl/©m^^KM-ra^i«l©rS | Z(m)- 
Y, Cm) II 1 i II Z(m) + Y, (m) I '4it#lt. (5)^©*IJS 

fkttimtcim* y, *s*£ft&fc&sft&i> c i^*^ 
a. 

[0026] (5)S-r©»ff«. (2>^^fijffl-ra =fcO 

4. ft»«©.*.-C«!)*WT*«. £IS. (2)^fffiBT 
*fca©tf#«4«*0T**£» aa^l^tt. 40 

hJU©^tt^©»»iJPllWifiW«:tta. -#(5)S 

©ffffitcB. ap»anB©^7aR»©«#ifti*"c+^r 
wnB/jNStHJi'. n-m«©6»*iS8^ai*a. 

[ 0 0 2 7 ] «T. 0ffl*ffl^-C*»W© 1 

SBfflKSMB-f*. 

[0 028] (JUfeflSSS 1301 »***»•©"?*- * 

«*SB©«Wi*jS'r. Amaian o n** afnajg 

^+-#-h-^©^7 ; -i»-f'3v>K?:A^-ra^S 50 



★ S(Z,Y,)^ lY, II 

mmtc 

5(-z,y,)^ By, I 2 - II pcz+y,) V (4) 



I Y, II 2 - II P(Z+Y, ) 6 2 ) (5) 

fiSO . F-a^A^sftfcA^-* *rtain-» 
i o 2 ^iift u . rtsimss io2*^ 6£fi$ n 

[0 02 9 ] 3fflllMtc$tt£hTl>*afflttft; 

0 2fCjt{fU > SHKlStt«KiKl 0 2*>6£ffUc-f- 
*4ilfilli««:^l/"C«Sl$*ifcT F uxtcajfi-r*. 

[003 0] <s(Kett»$ia l 0 2 Sffi^H- 

^«gl 0 2©^S^Ii?r@2©7 0-^ + -F?rffll> 



(6) 



[ 0 0 3 1 ] ST. Xf ? 7"S20i-CAUi*^S 10 1** 

Ol*m<0&L. fcfttfXf *:7S2(B^iitf. Zt-j?'S2 

02T«, A#f-2#T-£^~X©f-££M$?rTS 
*>©*>£" 5 -e5T2>ft«xf 7 7"S203^ 

* 5 -eft if n« x f 7 7's204^iity. x f 7 ys203r 
«. sa-raiiM!Bi**tTU -ena^TUc^xf 7 

7'S201^ jttf . Xf 7'S204f «f£A;>3 f - 2 

S© *>©*> E *> b . -5" 5 rfcfttfX f 7 7*S205 

-v, -e-^T&WtltfXf ? T^Ol^jttf. Xf y^OS 10 

f 7 7*S20i-^iitf . 

[0 0 3 2 ] 094ffl<,>tXf- y7 , S2O3'C||tf?n^mf 

[0033] imss-ck. 7'a f 5? ^ :/jb£kht s ? 
a h ^-Y7"©^Ji'A(D2St. #J&5Hi(cM3~£v-f 

?a v* 4 ■fxm^zffltobtc 0 x k it,> 5 2 as© 

[0 0 34] Xf7 7'S80rCn£ 1 (CS^f £. s> 

{ (Y. (n) , n) , Of, (n) , n) , • • • , (Y. (n) , N) } 

Xf 7 7'S803r±ia^75rn^^fii©^)itC^#X. 

So 

[0 0 3 5 ] {(YIn(l)(n),In(l)),(YIn(2)(n),In(2)), 
••• ,(YIn(N)(n),In(N))} 30 
fcTF. ^©©StftlDg^&tfOJ^VX hVh. 
Ini-Ti : 

V„ = {Y,„ (1) (n),Y t . (1) (n), ••• ,Y,„„(n)} (6) 
I„ = ft (I) .1.(2). ••• .1.00} (7) 
Xr f7'S804r«n?r^>^U^> h 0. F-Jl^X 

v vsmi^m. 2as© x h ©w^»^©i fc o t 

* & . V, ©5teSK#> 6 n+l# @ ©fflV, (n+l) « n#g ©fitv 

, (n)feLht?*-5. *fc. I f (n)5rI(Sf iTST'o h £ 
-Y 7'YI, (n)©, JfcftfilYI, (n) (j)#V, (n) t — gtt" Z> . 40 
[0036] ^9.ttB3*fflC»TXf-»:/S20yt?»TS 

tih&mmmzmitz. i/e. 

<Dtc£>(D'<i> \-)lX (feTFf X F-^7 h^iPf^) 
[ 0 0 3 7 ] X f v 7"S30lCttf X h < £ h Mtltfif 

Z=GX 

xf 7 7*s302-C«-<^ MU2H©SJg©4 >f 2X 'J x 

vzmtz. ctummmmmz&mtzmm so 



2002-318818 
©Ji#£S©£ 'J X F- r& 0 . flRtf. rXh-c^hJl/ 

x©/$*Mi©*&*Kii©A & z> yxhT* 

£: 

A={A, ,A, .•••.A d } . (8) 

*/c/h$CH3:^^6Lfi©S^ ; &A l <!:S< : 

A t = {A 4 .i,, .Aj.l, ,A.} (9) 

HKmfcAiittaS-fS. 

[ 0 0 3 8 ] X f 9 7"S303r»FTR&t>'IWj!tT &£$© 

vimwmznmh. cfDsmzmtzm^-cmiT 

So 

[0039] Xf 7 ^S40irmfiS7>ffi© 7-f -f > ^ • 

•;x F-v.£flW#1" £. xf 7 7*s402r«f-x v*<t> 

©m^TKIX^CCSfejai^fil^V.^^^l/, *©fi[g 
^PTR-tCtSSM-rS. EP^ : 

IV. (PTR' )X(m) | ^ IV. (i)X(m) | , V, € {1,2 , • • • ,N} 
PHitc. -X(m)fC®^a^ffl^V„^6^U. ^©tta 

|V. (PIR- )+X(m) | ^ |Vm(i)+X(m) | , V jE {1,2, • • • ,N} 
Xf 9 7-s403-C«. BBii^4SBRti«T©J: ^KfiJJWft 
ft, : 

FTR* t = FTR* , ,BND' U =0,CAL\ = 0 
FTR\ = PTR\ , ,BND\ = 0,CAL'„-0 
PTR- l =PTR- 11 BICT l =0,CAL\=0 
FTR" „ = PTir , BNU „ = 0 ,CAL" „ = 0 
FTR= FTR* 

xf7 7's304f ttmvftmm&oxnmtefi ^. c© 
temzm 5 sfflorttwr s. 

[ 0 0 4 0 ] X f » ^"SSOl-CttjfifiNR^N, (X)*S*'& 

f 7 7'S503T«. 06 £fflt>T&jfrf -SFTROfgilfSrtT 

^ 0 Xf 7 7'S504r«IDS : #I. (PTO)©"X'D h *> ■{ 7\ 

..cp,,, iA.^rlt^T -S : 

S= I. (FTR) 

A, = 2 |p(Z,Y)l 

M6b±^0^m i att^ffi-en-esnri^© 
-c. gett«a^6©s*ttJi//c^-ca^. 

[ 0 0 4 1 ] X f -7 7"S505T«a^^N t .,001c ID§ 
^iA.©ffli*iMi-rs : 
N t (X) = N t .,()0 + {(s,A,)} 

xf •7-x , s506f tzjyfi'sy Ft. k*iin«xf 

•7 7'S507^m&. Z *> -CfeW*a«Xf 7 7'S503^^S. 
[0042] X f 7 7*5507-C«i5<5l*^©^^'eA, © 
*/JMi* Si-i. S/c^n(C>Ptfc;-r SlD§#s£ T t .i i 
LtEtt-Ti. JiLh-CXf 7 7*S304©jfi®ll£©SMt 

[ 0 0 4 3 ] Xf -7 7'S305T«t*k+l«:S5£1' Z>, xf 
7 7*S306T ttPTRCHSf^Htf t S. C ©MS^rH 6 . 

I 0 0 4 4 ] Xf -7 7*S60rCttBNU t ^OPTRT^S 
O^i'^^^20. ^.t*ntfXf7 7"S602^ ^ 



' • • (7) 
11 

^r&WfttfXr 7 7"S603~-iIt?. 7'S602t , «U( 
T©A&a£flirtS : 
BND\ = l,Dx~ L =°° 

X r v 7'S603T BND" , ft 0 *>0 PTR" „ #PTR* „ £: 3? U I » 

fttlftZtf&U. ■£5?*fttfXf-y:7'S604'x -e^ 

7'S605^jlty„ X^7 7'S604-C«Ji(T 
BND* n = l,Dx-„=<» 

Xf7 7'S605T«BND* L A* 0 #>OPTR* L ^FTR-Hct 0 /h§ 

<,>£'ofr£$2U •£5-c<fcfttfXf--7-/S606^ t -e^ io 
-c&tthtfxf- 7 7s607'Njltf. 7\f~v y's606Vf,tUTf 

BND" „ = 1 , Dx~ „ = 00 
BND' L - 1 , Dx' t = 

X r 7 ?'S60TV « BND* „ # 0 #>0 PTR' „ #N2 W±#> <t* *5 

t*x f - 7 7s609^jttf „ x f - ? -/seos-cteWT©^* 

fSfTTS : 

BND* H =l,Dx*,,=°° 

Xf 7 7'S609T«BND- L <b BNO" „RO'BND\ £ BND\ ©4 11 20 

^71, ^r&WfttfXT^T'seiO'Njltfo 

[0 045] Xf- ^7 , S610T«BND- L +CAL\*H #>£'5 

Dx\=|V.(PTR- L )-Z(m)| 
CAL\ =1 

X f - 7 ^Se^-Cli BND" „ +CAL" n ftlft£? #>£&2 1/ . 
-e^T*n«Xf v-?'S613^. ^^WniiXf^ 30 
S614^jlif„ Xf"? 7'S613rttWT©^*^T"$"S : 
Dx |V.(PTR- M )-Z(m)| 
CAL',=1 

7,7- -J 7' S614T (J BND* L +CAL' L ft 1 ft £ '*> ft £tfc2 1/ » 

-e- ^ -c^ntix f » 7*s6i5^. -e 5 -e&tttitfx r » 7' 

SSifi^iitJ. Xr-V7*S6l5T*«WT©®a?:^fTS : 

Dx\ = ZV„(PTR' L )-Z(m)| 
CAL*i =1 

7'S616T I J BND* , +CAL* „ ft I ft £' 9 ft £t&2 1/ . 

%'>-C8>i\&ZT'j 7"S6i7^, -e^T^t-rtvKXf- -7 7* 40 
S6i8^jitJ„ ^-f »^S6im*felT©«Hl*lltf"rS : 

Dx*„= |V„(PTR'„)-Z(m)| 
CAL'„=1 

XT-^SeiBTliDx-^Dx-^ Dx' L . Dx\©(,>-f tufc 9 
fc/M Wfltf * f 1 v 7*S619^. * 5 "Cfctftltf* TV? 
S620^jt&. 

[ o o 4 6 ) x f -j y'ssi^mromm^my u x 

r v 7 , S306©PTR©i!$r : £ifi*7"3~5 : 
Dx= Dx" L , PTR= PTR" t ,CAL" L = 0 

X^^7*S620-C»Q>r»^DX- v . Dx* L> Dx*„©^-fft<fc 0 50 



^2 002-3 1 8 8 1 8 

.12 

4>/Jn$ Wfttf Xr 7 7*S62l^. ^ orfcttntf Xx -j V 
S622^Mts. 

[ 0 0 4 7 ] X f 9 7S62lTttttT©*Ml*£ff O. X 
?--7 7S306©PTR©H&r£il&7'f 5 : 
Dx=Dx" B ,PTR=PTR" M ,CAL" h = 0 

ZT vy'S622VitDx\ftDx\, Dx"„, Dx*,,©t>-fflJ: <3 

fc/M*wft'Jxf--7 7 , s623^. *9-e&i**i«xy » 7 - 

S624^jttf. 

[0048]Xf-; 7"S623'T«WT©toa«:JlTf U X 
f - » ^oeCPTR©!!^?:*^^ £ : 
Dx=Dx\ ,PTR=PTR* 1 . ,CAL\ = 0 

Xf'j 7'S624"C«Ji(T©^a?: Htf b . X t -j 7*S306© 
^©MSf^TI" -5 : 
Dx=Dx*„ ,PTR=PTR*„ ,CAL*„=0 
PTR©Mirtea-C«. PTRiMil-r4*»O«[*gCM0. 

»7*T4. 

[0 049] Xr^7'S307r«, ^S-Cp*Sf>'p-?:it 

#-rs : 

P*= lY. l-IZl 
P" = I Y, 0 - 1 z 0 

r (s) 

[00 50] Zf-v 7*S309*C«, Xf 7 7*S308rtf o7c 
Ma©K ») fit* 5 FALSER i' yftZnmb, FALSER 6 X f 
v 7 - S310^ . * ^ V & X 7 1 - 7 -/S312^3ltf . 
[ 0 0 5 1 ] Xf-7 7"S310-C«, &®m£ftt> € t -,tC 

jp-rs : 

N t (X)<-N t . 1 (X){(T t .,,^ t . ,)}+{(I. (PTR),Dx)} 
Xr 7 7*S311-e«. ^(^©^^©^^©DxWS/hti?: 
|,{C, *©ID*#*Tc«:t&ttU Xf-; 7'S313^it 
tf. 

[ 0 0 5 2 ] X r 7 ^ss^-Ctt^S^H^ L/, X f 7 V 
S313^jtt? : 

Xt" 7 7'S313t?«t* 'f>i'b7<>ht> 

mz.tct>mmzm * vt&wntfxf? 7's306^ 

It?. 

[0053] Xf 5 7'S308©|lJS^a'CllTf-r€>H^lt 

[ 0 0 5 4 ] Xt" 7 7*s7oir«.i?ri(cg:5S-r 2>. Xf" 7 
7'S702-C«WT©^a*Hif-r?> : 
p'=p'-(Ys(A,)-X(A 1 ))* 
p^p'-CYsCAO + XCA,))' 
Dx'-Dx-CYsCAO-XCA,)) 1 

x f 7 7s703-m:X5££¥iJS O . ffiS U It hm 0 ffitcTR 

iE*a9SEo-c«ia*»7o. * i 7-cA c ct-tn«xf-7 7-s7 

04^3tt?. Xf 7 7"S704r«. j4 -f > ^ > h U 



13 

^©asu^? Y>\>o'#jtd&m7Ltch. mvmc false 

[005 5] T.rv 7'S3l3-Ci^7 Ucif©, N, (X)* 5 

[oo56] uiwnucmmm^c^^^nnm, 
wmtc£^x®w?z>. 

[0057] am^mm ±&LtcmmBm<Dm)& 

Y 7"?fck= 1<M. 7'd h £ ^ :/»n= lOOOOUKft 

• Y)\><D'#jt : d={10#>P>100£"C\ 10#£\&} 

• CPU : Penti umlll , ( 500Mte) 

• y -f>^t :128MB 

• OS:Linux-2.0.36 

ft, 7'd ^7 Affile zm^tc. 

[0 05 8] CH^HD 

(2) -m^ffl^T^ d&7C ^2 hA2IHI©tW? , > 

(3) -«S**JBl>T. <«X7C©f-X h^* MU* 1 
C4) ±S5[^^^tT^. 

[0 05 9] ±§e5o©3Mi£ioJn|iSii)jIU 
ffl*fcpuu#ra= (ff&T Jb=f 'J XA©CPlJ5m)/C&R8£ 

©crjsih) 

HH^S^l 1 0 fC^T. hJKD^7C(numbe 
r of dimensions)^!! tl^CPUB#l9 (CPU time r 
atio)««tt» y'uYm^mZ'^JZtbXmb 
tc. 

[0 06 0] @J;97"Ph£0'ifcNfCMfc>6-f, ^9 

h wATvommicw -> x mt c p u usra # 1 - 
2"- rff UP i> s c t tft>fr i> . i/A i L/N= iooooH© 

£ * ©ffi * tt#&tC /J> 3 ^. C © t £ CPUB#Ptg© 
ffif», 10^7C"C2% lOO^TCffelS^i^ i^#7cTfe+# 

[006 1] <I»;<§2 >±ai/fcS«H«««*B» 

[ 0 0 6 2 ] m 1 1 «. *H*»J»©BS^K©flteS% 

^r-^$n/cNfi©^a h # ^ ^©n^isttsnT 

[0 06 3] f-^-Xl 0 3(C«iWSn^d^7C^ 



(8) #P8200 2-3 1 88 1 8 

14 . 

A,. •••) ©m$©£iK. SitiHWcfcWSBiflMc. *© 
A«W>tt«*H«OTtf£C£KJ:^T. *JdE-T4Ci3W 

A^BO«ijfi^kfli©ii«*ijajR$n. -en-enccng^ 

tfwe>ft/tAfe#-st©i«K£J6ftU «K©*feJSi»A 
»*. H&K»J£Tr£A^T*££^£T£c<t#ffi* 

[0 0 6 4 ] *fc. SBMSttgLgg 1 0 2©MS^« 
io ». m&Bt&l £|§M§cc02«:Se5 o {bu XT-j^Si 

04~205-C«. Mtt, ^&J!ll>fc38ia&9£% 
[ 0 0 6 5 ] m 1 2 CCfcOT, Xf J 7'S301~313«, 

03<tp]-r*6. ^7"^y's3i4Vkt, ftwrnsicmt 

4. X^yyS315-C«> d»7*##©*(K©**l f »IB 

fect^u. «K«{cBfjeffl**t»»^au©i'7^*i8» 

[o o e e j fisa. rntmt. *-©«S*»6a*i8IK: 

itfflb-cfe, *»©«ws*6««3ti4s/^^Atcaffl 

30 ctl>. 

[0 06 7] HOC, gg*^tH*'>^'?Art©3>h , > 

©^S(Cct0, ±ifi©«l6*KP13n4»^fc«*n 
[0 06 8] Cft6©tt£\ Wa^An- F^IB 
40 [006 9] 

h;b©rta©^fiitca-^<^ma ; 5r> isatcn-jf 

-T4C4*«t?*e©t?. C©SHHaKX-3l»T. 
^^^H -T 4 9 7 ^ ©BSJt*«a«:^tf 1" £ C <h *i"C * 

naffi©jfflm^siw] 

[01 ] *#BB«:«4-^«9^»©^-f1ft^B©* 
so [02] MfHWUi©?M«5%-*- ^ a - ? + - h xh 



(9) 



15 



So 

[05] ^mm-^oymmit^m^t 7a-?t-n 

[06] Pra©ISIr****"r7D-^*- h-c*s. * 



12002-3 1 8 8 1 8 
16 . 

[07] mtDm&^MiTntJB-?-*- hV$>Z Q 

[08] fij^a^ii^^-r^cj-?-*- h-c&&„ 

[09] HUiffiKD^Ii^T^O-? + - 

[010] 9mmma*A*w*avii>*. 

[011] BBI»B©<«RJE**r^o » *0T?*S. 
[012] Stt*affl©»MI#«*mT70-^ + -hT 

[013] mrbJ yxAfcBMM-sHr**. 
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